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EXPERIMENTAL  RANDOM  FATIGUE  IN  ELASTIC  RANGE-FIRST  ORDER  MODELS 


T.  C.  Huang  and  Vinod  K.  Nagpal 
Department  of  Engineering  Mechanics 
University  of  Wisconsin-Madison 
Madison,  Wisconsin  53706 


Abstract 

An  experiment  program,  based  on  probabilistic  parameters  and  experi¬ 
ment  design,  was  conducted  for  the  fatigue  life  under  random  vibrations. 

A  batch  of  24  specimens  were  tested  to  study  the  effects  of  8  proba¬ 
bilistic  parameters  on  the  fatigue  life  in  3  designs.  The  response  of 
each  tested  specimen  was  fitted  with  time  series  models  from  which  the 
spectrum  and  spectral  moments  were  computed.  The  magnitude  levels  of 
all  8  probabilistic  parameters  for  every  specimen  were  computed  from  its 
spectral  moments  and  were  coded.  The  first  order  models  to  predict 
fatigue  life  were  obtained  by  regressing  the  log  of  fatigue  life  on  coded 
probabilistic  parameters.  From  the  table  of  analysis  of  variance  the 
F-ratio  is  computed  to  check  if  the  model  is  acceptable.  The  tables  of 
predicted  lives  together  with  residuals  and  95%  confidence  intervals 
are  also  given  for  each  model.  The  percent  deviation  from  the  actual 
life  of  the  predicted  life  obtained  from  the  first  order  models  of  life  _____ 
predicting  equations  were  in  general  comparable  with  or  better  than  the  1  j 
conventional  deterministic  fatigue  tests.  Furthermore,  this  method  is  ] 

in  contrast  with  the  existing  method  which  is  based  on  linear  cumulative 
damage  and  cycle  counting  and  normally  involves  several  hundred  percent  ' 


error.  Therefore  the  existing  method  remains  academic  as  it  is  not  C0(j( 
suitable  for  application.  |fil3t 
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INTRODUCTION 

The  principle  of  linear  damage  accumulation  based  on  equilvalent 
cycle  counting  has  been  repeatedly  proved  to  be  unreliable  to  predict 
fatigue  life  under  random  loading  [1-12].  These  references  [1-12]  are  a 
selected  group  which  have  been  discussed  in  [13].  In  applying  this 
linear  cumulative  damage  principle,  the  random  loading  is  evaluated  in 
terms  of  equivalent  deterministic  cycles.  The  equivalent  cycles  for  a 
random  loading  have  been  obtained  by  several  investigators  using  differ¬ 
ent  methods  of  cycle  counting  [13,14].  Each  equivalent  cycle  contributes 
its  share  to  damage  based  on  deterministic  tests  of  that  amplitude.  In 
the  principle  of  linear  damage  accumulation,  the  shares  of  damage  due 
to  all  the  equivalent  cycles  are  sumed  up.  The  fatigue  failure  is  pre¬ 
dicted  when  the  damage  sum  reaches  a  value  of  1.0.  The  life  is  con¬ 
sidered  to  be  overestimated  if  the  actual  sum  remains  below  1.0  and  is 
underestimated  if  the  actual  sum  exceeds  1.0  at  failure. 

The  above  mentioned  representative  references  dealing  with  the 
linear  damage  accumulation  principle  have  reported  deviations  as  follows. 
The  factor  of  overestimation  of  the  fatigue  life  using  the  peak  counting 
method  ranged  from  1.5  to  3.5  in  [1],  from  2.0  to  3.5  in  [2],  and 
reached  a  value  of  10.0  in  [3].  Kowalewski  [4],  Locock  and  Williams  [5], 
and  Brown  and  Ikegemi  [6]  have  used  the  mean  crossing  peak  count  method 
and  have  found  that  the  fatigue  life  was  overestimated  by  a  factor  which 
went  up  to  3.0  in  [4],  ranged  from  -5.0  to  15.0  in  [5],  and  varied 
between  2.0  to  3.0  in  [6].  Dowling  [7]  used  the  rainflow  counting  method 
and  reported  that  the  predicted  fatigue  life  in  some  cases  was  as  high  as 
3  times  the  actual  life.  The  following  investigators  have  used  counting 
methods  to  evaluate  the  fatigue  loading  in  equivalent  cycles  but  the 
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methods  were  neither  mentioned  nor  described.  Schultz  [8]  reported 
that  the  fatigue  life  was  overestimated  by  a  factor  up  to  10. 0,  and 
Buch  [9]  found  that  the  life  was  overestimated  by  a  factor  ranging 
between  -2.95  to  6.0.  In  Smith  and  Malme  [10]  and  Swanson  [11]  it  was 
mentioned  that  the  fatigue  lives  were  overestimated  but  the  factors  of 
overestimation  were  not  given.  On  the  basis  of  their  study  Ford,  et.  al. 
[12]  mentioned  that  the  failure  that  most  frequently  occurs  under  a 
general  loading  spectrum  is  not  predictable  from  the  constant  amplitude 
tests.  The  counting  methods  and  the  factors  of  overestimation  obtained 
by  several  investigators  using  various  counting  methods  are  summarized 
in  tabular  form  in  [13]. 

The  review  above  shows  that  the  principle  of  linear  damage  accumula¬ 
tion  is  not  accurate  to  predict  the  fatigue  life  under  random  loading 
and  is  therefore  not  reliable  for  industrial  use.  Recently  a  new 
history-dependent  stochastic  model  of  cumulative  damage  is  being 
developed  by  Bogdanoff  [15,16,17]  by  taking  a  comprehensive  view  of 
the  entire  failure  process.  The  model  includes  major  sources  of  vari¬ 
ability  encountered  in  the  complex  fatigue  process  by  acquiring  data  of 
a  specific  type  in  addition  to  the  data  on  time  to  failure.  This  data 
is  essential  for  improvement  in  predictive  accuracy.  It  also  implies 
that  the  frequent  inaccuracies  encountered  in  life  predictions  when 
using  current  cumulative  damage  models  cannot  be  substantially  reduced 
using  such  models.  In  this  report  an  entirely  different  phenomenological 
approach  has  been  undertaken  and  an  experiment  program  based  on  8  proba- 
listic  parameters  was  conducted  using  experiment  design  technique  [18]. 
The  models  were  developed  which  would  produce  a  reliable  estimate  of 
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fatigue  life  for  the  given  levels  of  the  probabilistic  parameters. 

Since  the  purpose  of  this  research  is  to  develop  a  novel  methodology 
for  random  fatigue  based  on  the  probabilistic  parameters  and  experiment 
design,  the  choice  of  the  specimen  materials,  types  of  loadings  (axial, 
shear  or  bending),  and  range  of  parameters  is  irrelevant.  Therefore, 
they  are  arbitrarily  chosen. 

A  comprehensive  collection  of  references  related  to  random  fatigue 
prior  to  1968  was  provided  by  Swanson  [19).  Bogdanoff  [15,16,17]  has 
highlighted  many  references  particularly  in  the  probabilistic  and 
statistical  aspects  of  cumulative  damage  in  developing  his  new  cumula¬ 
tive  damage  model. 

I.  EXPERIMENTS,  PARAMETERS,  DESIGNS  AND  MODELS 

1.  Experiment 

The  experiment  program  was  conducted  with  24  specimens  numbered 
from  1-24  in  3  designs.  A  bending  specimen  of  aluminum  alloy  6061-T6 
was  chosen  for  the  experiment.  The  same  specimen  and  the  experimental 
setup  were  used  in  a  preliminary  investigation  of  fatigue  failure  of 
materials  under  narrow  band  random  vibrations  [20,21].  In  [20]  the 
specimen  and  the  experimental  setup  are  described  and  pictorially 
illustrated.  The  dimensions  of  the  specimen  are  so  designed  that  a 
constant  strain  is  obtained  all  along  the  test  length  under  any 
bending  load.  The  dimensions  of  the  testing  section  of  the  specimen 
are  shown  in  Fig.  1. 

2.  Parameters 

Fig.  2  shows  a  typical  random  response  signal  in  elastic  range 


which  is  characterized  by  8  probabilistic  parameters.  They  are  mean. 


variance,  zero  upcrossings,  level  upcrossings,  duration  of  excursion 
above  zero  level,  duration  of  excursion  above  level,  band  width, 

and  average  amplitude  above  level.  In  these  variables  is 

the  strain  level  corresponding  to  a  material  life  of  107  cycles  in 
the  deterministic  fatigue  tests,  and  is  the  yield  point  strain  of 

the  materials.  A  mathematical  analysis  of  these  8  probabilistic 
parameters  is  given  in  [22].  Hereafter  these  probabilistic  parameters 
will  be  simply  referred  to  as  variables  as  well. 

3.  Designs 

The  first  design  of  the  experiment  consisted  of  10  tests  and 
formed  a  full  factorial  design  with  2  center  points.  The  second 
design,  a  central  composite  design  with  four  center  points,  was  formed 
by  adding  8  more  tests  to  the  first  design.  The  third  design  was 
formed  by  adding  to  the  second  design  random  replications  of  six  tests 
from  the  first  and  the  second  designs. 

a.  Range  The  first  design  was  based  on  the  three  different 
magnitude  levels  of  the  three  variables;  namely,  mean,  variance  and 
zero  upcrossings.  Two  more  levels  from  each  of  the  3  variables  were 
added  in  the  second  design  so  that  one  was  lower  and  the  other  one 
was  higher  than  the  3  previous  levels.  In  the  third  design,  6  tests 
were  replicated  without  introducing  new  levels  of  the  variables. 

Thus,  the  second  and  the  third  designs  had  5  levels  from  each  of  the 
above  3  variables.  The  lowest  and  the  highest  of  the  3  levels 
represent  the  upper  and  the  lower  bounds  of  the  range  of  variables. 
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The  upper  bounds  of  the  mean  and  variance  are  limited  such  that  the 

response  remains  within  the  elastic  range.  Otherwise  the  range  of  the 

variables  can  be  arbitrarily  chosen.  The  ranges  of  the  mean,  variance 

and  zero  upcrossings  based  on  5  levels  are  110-1800  micro  inches/inch, 

2 

733,373  -  2,250,935  (micro  inches/inch)  ,  and  7.0  -  13.8,  respectively. 

b.  Code  The  above  3  variables  which  were  controlled  in  the 
experiment  were  coded.  The  magnitude  levels  of  the  remaining  5  vari¬ 
ables  which  could  not  be  controlled  in  the  experiment  were  computed  from 
the  response  of  each  test  specimen.  The  response  of  each  test  specimen 
was  digitized  and  recorded.  The  time  series  models  were  fitted  to  the 
recorded  response  of  every  test  specimen  and  the  spectrum  and  spectral 
moments  were  computed  from  the  time  series  model  parameters  [23].  The 
actual  magnitude  levels  of  the  variables,  were  computed  and  their  magni¬ 
tude  levels  were  coded.  The  variables  and  the  actual  levels,  and  coded 
variables  and  coded  levels  are  shown  in  Table  1.  Tables  2  and  3  show 
the  actual  levels  and  the  corresponding  coded  levels,  respectively,  of 
all  variables  for  all  24  tests. 

4.  Models 

The  first-order  regression  models  were  obtained  for  all  3  designs 
by  regressing  the  log  of  actual  fatigue  lives,  y,  of  the  specimens  on 
corresponding  coded  levels  of  variables.  The  regression  models,  which 
in  our  case  represent  life-predicting  equations,  predict  log  of  life, 
y,  of  the  materials  for  any  coded  levels  of  variables.  The  predicted 

A 

fatigue  life,  T,  was  computed  by  taking  antilog  of  the  y. 
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Two  first  order,  life  predicting  equations  were  obtained  for  each 
of  the  3  designs.  The  first  was  obtained  by  regressing  y  on  3  coded 
variables,  and  the  second  one  was  obtained  by  regressing  all  8  coded 
variables.  For  each  equation  the  tables  of  analysis  of  variance,  and 
of  predicted  lives  together  with  residuals  and  95%  confidence  intervals, 
were  constructed.  In  the  analysis  of  variance  the  F-ratio  was  computed 
for  every  life  predicting  equation  in  order  to  determine  whether  the  regression 
was  effective  and  the  model  was  acceptable.  The  confidence  intervals 
were  computed  using  the  standard  deviations  of  y  and  the  t  value 
from  the  t-table  with  degrees  of  freedom  equal  to  that  of  residuals. 

The  distribution  of  residuals  for  every  equation  were  carefully  studied 
for  any  trend  or  pattern  present  in  the  residuals.  In  case  of  any  trend 
or  pattern  in  the  residuals  the  model  was  considered  inadequate. 

II.  FULL  FACTORIAL  DESIGN  WITH  2  CENTER  POINTS 

The  two  life  predicting  equations  are  obtained,  using  first  the 
ten  tests  of  the  experiment,  by  regressing  the  log  of  life,  y,  on  the 
coded  variables  under  consideration.  These  ten  tests  form  a  full 
factorial  design  with  2  center  points.  The  first  equation  consists  of 
three  variables  which  are  mean,  variance,  and  zero  upcrossings.  The 
second  equation  consists  of  all  8  variables.  The  tables  of  analysis 
of  variance,  and  of  predicted  lives  with  95%  confidence  interval  are 
constructed  for  each  equation. 

1.  Three  Variables 

The  life  predicting  equation  is  obtained  as 

y  -  5.94  -  0.224x1  -  0.690x2  -  0.080x3  Cl) 

The  analysis  of  variance  is  given  in  Table  4.  The  F-ratio  of  21.65 


with  3  and  6  degrees  of  freedom  is  obtained.  The  corresponding  F  value 
from  the  F-table  at  95%  significance  level  is  4.76,  which  shows  that 
the  regression  is  effective  and  that  the  model  is  acceptable.  The 
residual  sum  of  squares  is  0.4897  in  comparison  to  a  total  of  5.7888, 
a  8.5%.  The  other  92.5%  of  the  total  sum  of  square  is  due  to  the 
regression.  Analysis  of  variance  indicates  that  the  zero  upcrossings 
ha\e  no  significant  effect  on  the  fatigue  life. 

The  predicted  lives,  together  with  residuals  and  95%  confidence 
intervals,  are  given  in  Table  5.  A  few  of  the  confidence  intervals  are 
quite  wide  because  of  the  small  number  of  tests.  The  residuals  appear 
to  have  a  sinusoidal  pattern.  This  indicates  that  some  variables  which 
may  have  significant  effect  on  fatigue  life  are  missing  from  the  equa¬ 
tion. 

2.  All  8  Variables 

The  life-predicting  equation  is  obtained  as 

y  *  5.90  -  0.226x^  -  0.687x2  +  0.519xg  +  0.261x^ 

+  0.378x5  -  0.382x6  -  0.285x?  -  0.509xg  (2) 

The  analysis  of  variance  is  given  in  Table  6.  The  F-ratio  for  this 
equation  is  4.796  with  8  and  1  degree  of  freedom.  The  corresponding 
F-value  from  the  F-table  at  95%  significance  level  is  239.  The  com¬ 
parison  of  the  two  F-values  indicates  that  the  regression  above  is 
not  effective,  even  though  97.6%  of  the  total  sum  of  squares  is  due  to 
regression.  This  is  because  the  number  of  variables  is  too  large  in 
comparison  to  the  number  of  tests. 

The  predicted  lives  (together  with  residuals  and  95%  confidence 
intervals, are  given  in  Table  7,  Most  of  the  confidence  intervals  are 
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fairly  wide  because  Che  value  of  t  used  in  computing  the  intervals  is 
very  high  for  1  degree  of  freedom.  The  sinusoidal  pattern  in  residuals 
has  disappeared  here. 

III.  CENTRAL  COMPOSITE  DESIGN  WITH  4  CENTER  POINTS 

The  full  factorial  design  was  augmented  by  8  more  test  runs  to  make 
it  a  central  composite  design  with  4  center  points.  This  design  has  a 
total  of  18  tests.  The  two  life  predicting  equations  obtained  using 
18  tests  are  given  below. 

1.  Three  Variables 

The  life  predicting  equation  is  obtained  as 

y  *  5.9008  -  0.1117x1  -  0.7291x2  -  0.0124x3 
The  analysis  of  variance  of  the  equation  is  given  in  Table  8.  The  F- 
ratio  with  3  and  14  degrees  of  freedom  was  found  to  be  33.86.  The 
corresponding  F-value  from  the  F-table  of  95%  significance  level  is 
3.34.  The  comparison  of  two  F-values  indicates  that  the  regression  is 
effective  and  that  the  model  is  acceptable.  In  this  case  the  effect  of 
zero  upcrossings  on  the  fatigue  life  is  also  observed  to  be  insignificant. 
The  residual  sum  of  squares  is  1.1841  in  comparison  to  a  total  of 
9.7843,  a  12.1%.  The  other  87.9%  of  the  total  sum  of  squares  is  due  to 
the  regression. 

The  predicted  lives, together  with  residuals  and  95%  confidence 
intervals,  are  given  in  Table  9.  The  actual  fatigue  lives  of  a  couple 
of  specimens  are  out  of  the  95%  confidence  interval.  This  indicates 
that  more  variables  are  needed  to  improve  the  prediction. 

2.  All  8  Variables 

The  life-predicting  equation  is  obtained  as 


9  =  5.823  +  0.039X,  -  0.9040xo  -  0.0817x_  +  0.1267x, 
12  3  4 

-  0.2310x5  +  O.lOOOxg  -  0.1118x7  +  0.0145xg  (4) 

The  analysis  of  variance  of  the  above  equation  is  given  in  Table  10. 

The  F-ratio  for  this  equation  is  computed  as  19.29  with  8  and  9  degrees 
of  freedom.  The  corresponding  value  of  F  from  the  F-table  at  95% 
significance  level  is  3.28.  The  comparison  of  the  two  F-values  indicates 
that  the  regression  is  effective  and  that  the  model  is  acceptable.  The 
residual  sum  of  squares  is  0.5394,  in  comparison  to  a  total  of  9.7843, 
a  5.5%.  The  other  94.5%  of  the  total  sum  of  squares  is  due  to  the 
regression.  The  analysis  of  variance  indicates  that  zero  upcrossings, 
duration  of  excursion  above  zero,  and  average  amplitude  above  level 

have  negligible  effectson  the  fatigue  life. 

The  predicted  lives  together  with  residuals  and  95%  confidence 
intervals  are  given  in  Table  11.  In  this  case  all  the  actual  lives  are 
within  the  predicted  confidence  interval.  The  confidence  intervals  are 
relatively  wider  in  this  case  in  comparison  to  the  case  of  3  variables 
because  the  number  of  variables  is  larger.  Consequently,  the  number 
of  degrees  of  freedom  associated  with  t  is  small,  which  results  in  a 
large  value  of  t. 

IV.  CENTRAL  COMPOSITE  DESIGN  WITH  4  CENTER  POINTS  AND  6  REPLICATIONS 

Six  more  tests  are  added  to  the  central  composite  design  by  randomly 
replicating  six  points  on  the  design.  All  test  numbers  and  test  points 
are  shown  in  Fig.  3.  Adding  these  six  runs  resulted  in  a  total  of 
24  tests.  Two  life  predicting  equations  were  obtained  for  this  design 
also.  These  equations  and  their  analyses  are  given  below: 
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1.  Three  Variables 

The  life-predicting  equation  is  obtained  as 

y  =  5.8545  -  0.1205Xl  -  0.7595x2  -  0.0063x3  (5) 

The  analysis  of  variance  of  this  equation  is  presented  in  Table  12.  The 
computed  F-ratio  is  52.48  with  3  and  20  degrees  of  freedom.  The  corre¬ 
sponding  table  value  of  F  from  the  F-table  at  95%  significance  level 
is  3.10,  and  a  comparison  of  the  two  F-values  indicates  that  the  regres¬ 
sion  is  effective  and  the  model  is  acceptable.  The  residual  sum  of 
squares  is  1.4781  in  comparison  to  a  total  of  13.1131,  a  11.3%.  The 
other  88.7%  of  the  total  sum  of  squares  is  due  to  the  regression.  Zero 
upcrossings  were  found  to  have  an  insignificant  effect  on  the  fatigue 
life. 

The  predicted  lives, together  with  residuals  and  95%  confidence 
interval, are  given  in  Table  13.  There  are  six  actual  fatigue  lives 
which  fall  out  of  the  predicted  confidence  intervals.  This  observation 
also  indicates  a  need  for  some  more  variables  to  be  introduced  to 
improve  prediction,  but  it  should  be  noticed  that  the  confidence 
intervals  are  relatively  narrow. 

2.  All  8  Variables 

The  life  predicting  equation  is  obtained  as 

y  «  5.826  -  0.037x,  -  0.9194x„  -  0.026 8x-  +  0.0782x. 

12  3  4 

-  0.1390x5  +  0.0530x6  -  0.0830x?  +  0.0437xg  (6) 

The  analysis  of  variance  of  this  equation  is  given  in  Table  14.  The 
computed  F-ratio  is  23.78  with  8  and  15  degrees  of  freedom  in  comparison 
to  F  value  of  2.64  from  the  F-table  corresponding  to  the  same  degrees 
of  freedom  at  95%  significance  level.  This  shows  that  the  regression 
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1.  Three  Variables 

The  life -predicting  equation  is  obtained  as 

y  *  5.8545  -  0.1205x1  -  0.7595x2  -  0.0063x3  (5) 

The  analysis  of  variance  of  this  equation  is  presented  in  Table  12.  The 
computed  F-ratio  is  52.48  with  3  and  20  degrees  of  freedom.  The  corre¬ 
sponding  table  value  of  F  from  the  F-table  at  95a  significance  level 
is  3.10,  and  a  comparison  of  the  two  F-values  indicates  that  the  regres¬ 
sion  is  effective  and  the  model  is  acceptable.  The  residual  sum  of 
squares  is  1.4781  in  comparison  to  a  total  of  13.1131,  a  11.3%.  The 
other  88.7%  of  the  total  sum  of  squares  is  due  to  the  regression.  Zero 
upcrossings  were  found  to  have  an  insignificant  effect  on  the  fatigue 
life. 

The  predicted  lives, together  with  residuals  and  95%  confidence 
interval, are  given  in  Table  13.  There  are  six  actual  fatigue  lives 
which  fall  out  of  the  predicted  confidence  intervals.  This  observation 
also  indicates  a  need  for  some  more  variables  to  be  introduced  to 
improve  prediction,  but  it  should  be  noticed  that  the  confidence 
intervals  are  relatively  narrow. 

2.  All  8  Variables 

The  life  predicting  equation  is  obtained  as 

y  -  5.826  -  0.037xx  -  0.9194x2  -  0.0268x3  +  0.0782x4 

-  0.139Qx5  +  0.0530x6  -  0.0830x7  +  0.0437xg  (6) 

The  analysis  of  variance  of  this  equation  is  given  in  Table  14.  The 
computed  F-ratio  is  23.78  with  8  and  15  degrees  of  freedom  in  comparison 
to  F  value  of  2.64  from  the  F-table  corresponding  to  the  same  degrees 
of  freedom  at  95%  significance  level.  This  shows  that  the  regression 
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is  effective  and  that  the  model  is  acceptable.  The  residual  sum  of 
squares  Is  0. 9585  in  comparison  to  a  Lota  I  ol  II. Mil,  a  I  .VI.  Tin- 
other  92.7%  of  the  total  sum  of  squares  is  due  to  the  regression.  In 
this  case  the  3  variables,  zero  upcrossings,  the  duration  of  excursion 
above  zero,  and  the  average  amplitude  above  were  also  found  to  have 

an  insignificant  effect  on  the  fatigue  life. 

The  predicted  lives, together  with  residuals  and  95%  confidence 
intervals,  are  given  in  Table  15.  It  is  seen  that  one  of  the  fatigue 
lives  is  out  of  the  confidence  interval  and  that  the  two  others  are 
on  the  border.  The  confidence  intervals  are  fairly  narrow.  This  is 
found  to  be  the  statistically  best  model  out  of  all  six  models. 

V.  DISCUSSIONS  AND  CONCLUSIONS 

An  experiment  program  based  on  8  probabilistic  parameters  and  experi¬ 
ment  design  has  been  performed  to  predict  fatigue  life  under  random 
loading.  The  regression  models  of  these  probabilistic  parameters  for 
all  3  designs  of  10,  18  and  24  tests  have  been  developed  to  predict 
fatigue  life  under  random  vibrations.  Unlike  the  linear  damage  accumula¬ 
tion  theory, the  approach  described  in  this  report  is  proved  to  be  reliable 
and  accurate. 

It  appears  from  the  analysis  of  all  the  models  that  the  following 
4  variables,  mean; variance;  duration  of  excursion  above  level, and 

the  band  width  showed  significant  effects  on  the  fatigue  life.  Among 
these  4  variables  the  variance  has  the  most  significant  effect, 
followed  by  band  width  and  mean.  The  duration  of  excursion  above 
ef  level  has  the  least.  The  regression  on  4  variables  reduced  the 
residual  sum  of  squares  by  5.6%  for  the  18  tests  and  3.6%  for  the 
24  tests  designs  in  comparison  to  the  one  obtained  by  regressing  on 
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variance  alone  as  shown  in  Table  16.  Any  reduction  in  the  residual  sum 
of  squares  results  in  a  corresponding  improvement  in  the  predicted  lives, 
y.  Consequently,  even  a  small  improvement  in  y  would  produce  a  large 
improvement  in  predicted  fatigue  life,  T,  as  the  models  predict  the  log 
of  fatigue  life,  y. 

A  comparison  with  respect  to  the  deviations  of  the  predicted  lives 
and  residual  sum  of  squares  of  all  the  models  investigated  in  this  report 
is  given  in  Table  17.  This  table  shows  that  the  model  of  all  8  proba¬ 
bilistic  parameters  gives  a  lower  residual  sum  of  squares  as  compared 
to  the  one  of  3  parameters  for  each  of  the  3  designs.  The  model  of  8 
probabilistic  parameters  for  the  first  design  of  10  tests  is  found 
inadequate  as  it  did  not  pass  the  F-test.  The  better  predicted  lives 
and  a  lower  residual  sum  of  squares  based  on  2A  tests  is  obtained  for 
model  of  equation  (6).  This  model  was  considered  to  be  the  best  fit 
model.  The  predicted  lives  of  the  best  fit  model  deviate  from  the  actual 
lives  within  a  range  of  -59.6%  to  24.1%,  with  an  average  deviations  of 
17.4%  on  the  negative  side  and  16.0  on  the  positive  side.  These  devia¬ 
tions  are  small  in  comparison  to  the  several  hundred  percent  obtained 
in  the  linear  damage  accumulation  theory. 

The  actual  lives  of  almost  all  the  tests  is  higher  than  the  lower 
level  of  the  95%  confidence  interval  of  the  best  fit  model.  So  the 
lower  level  can  be  considered  to  be  a  conservative  estimate  of  the 
fatigue  life.  The  predicted  life,  by  the  model  is  the  mean  expected 
fatigue  life  of  the  material  for  the  given  levels  of  the  variables. 

The  variables  which  showed  significant  effects  on  the  fatigue  life 
will  be  referred  to  as  significant  variables  and  considered  for  further 
studies.  The  analysis  of  a  few  more  first  order  and  some  second-order 


models  involving  significant  varibles  will  be  re-posted  in  [24]. 

VI.  SUMMARY 

Cl)  An  experiment  based  on  8  probabilistic  parameters  and  experi¬ 
ment  design  has  been  performed  to  predict  fatigue  life  under  random 
loading.  This  method  has  been  proved  to  be  a  reliable  and  accurate 
approach  in  contrast  to  the  linear  damage  accumulation  theory. 

(2)  The  experiment  was  conducted  in  3  designs  and  for  each  design 
2  first  order  life  predicting  equations  are  obtained.  The  tables  of 
analysis  of  variance,  and  of  the  predicted  lives  together  with  residuals 
and  95%  confidence  intervals,  are  constructed  for  each  equation. 

(3)  From  the  table  of  analysis  of  variance,  the  F-ratio  is  computed 
to  determine  whether  the  regression  is  effective  and  the  model  acceptable. 

(4)  Four  out  of  8  variables  have  been  found  to  have  significant 

effects  on  the  fatigue  life.  These  4  variables  are  the  mean;  variance; 
duration  of  excursion  above  level  and  the  band  width. 

(5)  For  each  of  the  3  designs,  the  models  of  all  8  probabilistic 
parameters  were  found  to  have  lower  residual  sum  of  squares  and  lower 
percent  deviations  in  comparison  to  the  models  of  3  parameters. 

(6)  The  model  of  8  probabilistic  parameters,  based  on  24  tests, 
was  considered  to  be  the  best  fit  model.  The  predicted  lives  obtained 
by  this  model  deviate  from  the  actual  lives  within  a  range  of  -59.6% 

to  24.1%, with  an  average  deviation  of  17.4%  on  the  negative  side  and 
16.0%  on  the  positive  side. 
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Table  1  Actual  Levels,  Coded  Levels  and  Coded  Variables 


UA 

„  ,  ,  .  Actual  Level 

Coded  Level 

Coded 

NO  i 

f di xauics  T 

Low 

Center 

High 

Low 

Center 

High 

van- 

ables 

l 

Mean  200.00 

447.21 

1000.0 

-1 

0 

+1 

X1 

2 

Variance  929500.00 

1284831.00 

1776000.0 

-1 

0 

+1 

x2 

3 

Zero 

Upcrossings  7.00 

8.97 

11.50 

-1 

0 

+1 

X3 

4 

Level 

Upcrossings  3.00 

4.24 

6.0 

-1 

0 

+1 

X4 

5  Duration  of 
Excursion 
Above  Zero 


Level 

0.65 

0.743 

0.85 

-1 

0 

+1 

X 

6 

Duration  of 

Excursion 

Above  ef 
Level 

0.20 

0.346 

0.60 

-1 

0 

+1 

x, 

7 

Band  Width 

0.935 

0.95 

0.965 

-1 

0 

+1 

X 

8 

Average 

Amplitude 

Above  ef 
Level 

2250.00 

2395.75 

2600.00 

-1 

0 

+1 

X, 
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Table  2  Test  Numbers  and  Actual  Levels  of  Variables 


Test 

No. 

Actual  Levels  of 

Variables 

1 

2 

3 

4 

5 

6 

7 

8 

1 

218.29 

880780.86 

6.94 

2.43 

0.57 

0.10 

0.9189 

2011.31 

2 

210.37 

930122.76 

14.00 

1.00 

0.56 

0.09 

0.9552 

2104.36 

3 

196.42 

1769478.37 

11.66 

6.91 

0.58 

0.18 

0.9730 

2386.04 

4 

1039.41 

934007.27 

7.30 

4.29 

0.85 

0.38 

0.9735 

2333.94 

5 

999.22 

1780656.42 

5.52 

4.01 

0.77 

0.75 

0.9192 

2712.85 

6 

1044.67 

1820240.94 

9.93 

7.29 

0.78 

0.40 

0.9674 

2729.46 

7 

963.97 

897936.15 

4.40 

2.37 

0.84 

0.31 

0.9300 

2313.13 

8 

210.12 

1792791.17 

7.07 

4.15 

0.57 

0.16 

0.9236 

2413.17 

9 

447.64 

1288513.58 

9.09 

4.01 

0.65 

0.18 

0.9698 

2183.95 

10 

436.46 

1274963.37 

8.92 

4.12 

0.63 

0.21 

0.9494 

2228.45 

11 

445.35 

1296591.74 

9.26 

4.65 

0.64 

0.20 

0.9548 

2259.02 

12 

465.25 

1299401.32 

9.39 

4.68 

0.65 

0.20 

0.9484 

2193.31 

13 

1800.23 

1255893.53 

2.81 

4.33 

0.94 

0.62 

0.9472 

3303.22 

14 

109.18 

1287049.27 

9.25 

4.39 

0.53 

0.13 

0.9559 

2283.18 

15 

460.22 

2371494.62 

9.14 

4.73 

0.61 

0.29 

0.9602 

2692.97 

16 

454.59 

716608.00 

8.63 

4.44 

0.67 

0.15 

0.9517 

1949.01 

17 

452.13 

1322376.67 

3.37 

6.31 

0.71 

0.11 

0.9556 

1973.89 

18 

453.28 

1394215.78 

13.16 

2.36 

0.63 

0.21 

0.9543 

1931.09 

19 

112.24 

1333860.83 

10.06 

6.37 

0.53 

0.15 

0.9859 

2193.51 

20 

202.90 

1768328.57 

11.40 

6.84 

0.58 

0.18 

0.9719 

2316.86 

21 

999.67 

1761848.86 

9.47 

7.16 

0.77 

0.39 

0.9772 

2814.93 

22 

199.47 

1763345.15 

6.99 

4.00 

0.56 

0.18 

0.9246 

2359.52 

23 

998.08 

950678.53 

7.42 

4.28 

0.85 

0.35 

0.9757 

2330.46 

24 

1008.71 

1755965.70 

5.27 

3.96 

0.76 

0.38 

0.9258 

2773.35 
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Table  3  Test  Number  and  Coded  Levels  of  Coded  Variables 


Test  _ _ Coded  Levels  of  Probabilistic  Parameters 


No. 

X1 

x2 

X3 

x4 

x5 

X6 

x7 

X8 

1 

-.891 

-1.166 

-1.035 

-1.609 

-1.979 

-2.355 

-2.100 

-2.551 

2 

.937 

-.998 

1.792 

-.161 

-2.111 

-2.454 

1.314 

1.486 

3 

1.022 

.989 

1.056 

1.409 

-1.849 

-1.192 

-1.136 

-.188 

4 

1.048 

-.985 

-.831 

.034 

1.000 

.168 

1.555 

-.493 

5 

0.999 

1.027 

-1.957 

-.162 

.234 

1.406 

-2.079 

1.588 

6 

1.054 

1.076 

.409 

1.562 

.388 

.262 

1.157 

1.672 

7 

.954 

-1.107 

-2.871 

-1.678 

.938 

-.202 

-1.339 

-.617 

8 

-.939 

1.029 

-.960 

-.066 

1.979 

-1.406 

-1.777 

.032 

9 

.001 

.009 

.053 

-.087 

-1.000 

-1.192 

-.032 

-1.412 

10 

-.031 

-.024 

-.023 

.263 

.867 

.911 

.327 

-.768 

11 

.005 

.028 

.127 

.283 

-1.116 

-.973 

.0988 

-.945 

12 

.049 

.035 

.183 

.057 

-.966 

-1.000 

-.179 

1.353 

13 

1.731 

-.070 

-4.677 

.098 

1.798 

1.060 

.400 

4.311 

14 

1.752 

.005 

.127 

.313 

-2.522 

-1.784 

.6842 

-.795 

15 

-.036 

1.893 

.075 

1.145 

-1.474 

-.032 

.380 

1.486 

16 

.020 

-1.803 

-.157 

-1.687 

-.774 

-1.524 

.293 

-2.986 

17 

.014 

.089 

-3.940 

-4.158 

-.342 

-2.088 

.354 

-2.811 

18 

.017 

.252 

1.543 

1.173 

.213 

-.885 

2.357 

-3.114 

19 

-1.718 

.116 

.461 

.584 

2.564 

-1.560 

.121 

-1.352 

20 

.982 

.987 

.965 

1.376 

1.894 

-1.222 

1.451 

.595 

21 

1.000 

.975 

.218 

1.511 

.292 

.202 

1.796 

2.099 

22 

-1.003 

.978 

-1.006 

.171 

-2.718 

-1.243 

-1.708 

-.343 

23 

.998 

-.930 

-.765 

.024 

.956 

-.007 

1.698 

-.514 

24 

1.010 

.965 

-2.144 

.199 

.166 

.158 

-1.626 

1.893 

Table  4.  Analysis  of  Variance  of  10  Tests 
First  Order  Model  of  3  Variables 
Life  predicting  equation: 


y  =  5.94 

-  0.224x1  - 

0.690x2  -  0. 

0600x3 

Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

F-Ratio 

Due  to  Mean 

0.3337 

1 

0.3337 

Due  to  Variance 

4.8717 

1 

4.8717 

Due  to  Zero  Upcrossings 

0.0937 

1 

0.0937 

Due  to  Regression 

5.2990 

3 

1.7663 

21.65 

Residuals 

0.4897 

6 

0.0816 

Total 

5.7888 

9 

F-ratio  is  greater  than  the  table  value  4.76  with  3  and  6  degrees  of 
freedom  at  95%  significance  level.  So  the  regression  is  effective 
and  the  model  is  accepted. 


Table  5.  Results  of  10  Tests,  First  Order  Model  of  3  Variables 
Life  predicting  equation: 


Table  6.  Analysis  of  Variance  of  10  Tests 
First  Order  Model  of  8  Variables 
Life  predicting  equation: 


y  =  5.90  - 

0. 226x1  - 

0.687x^  +  0. 

519x3  +  0. 

261x,  + 

4 

0. 378x 

5  -  0.382x6 

-  0.285x?  - 

0.5Q9Xg 

Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

F-Ratio 

Due  to  Mean 

0.334 

1 

0.334 

Due  to  Variance 

4.872 

1 

4.872 

Due  to  Zero  Upcrossings 
Due  to  Level 

0.094 

1 

0.044 

Upcrossings  Above 

Due  to  Duration  of 

0.203 

1 

0.203 

Excursion  Above  Zero 

Due  to  Duration  of 

0.006 

1 

0.006 

Excursion  Above  e 

0.073 

1 

0.073 

Due  to  Band  Width  f 

Due  to  Average  Amplitude 

0.057 

1 

0.057 

Above  £f 

0.005 

1 

0.005 

Due  to  Regression 

5.642 

8 

0.705 

4.796 

Residuals 

0.147 

1 

0.147 

Total 

5.789 

9 

F-ratio  is  smaller  than  the  table  value  239  with  8  and  1  degrees  of 
freedom  at  95%  significance  level.  So  regression  is  not  effective  and 
the  model  is  not  accepted. 


i 


Life  predicting  equation: 


24 


ri 


Table  8.  Analysis  of  Variance  of  18  Tests 
First  Order  Model  of  3  Variables 
Life  predicting  equation: 

y  =  5.9008  -  0.1117x1  -  0.7291x2  -  G.0124x3 


Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

F-Ratio 

Due  to  Mean 

0.1212 

1 

0.1212 

Due  to  Variance 

8.4720 

1 

8.4720 

Due  to  Zero  Upcrossings 

0.0070 

1 

0.0070 

Due  to  Regression 

8.6002 

3 

2.8667 

33.86 

Residual 

1.1841 

14 

0.0846 

Total 

9.7843 

17 

F-ratio  is  greater  than  the  table  value  3.34  with  3  and  14  degrees  of 
freedom  at  95%  significance  level.  So  the  regression  is  effective 
and  the  model  is  accepted. 


Table  9.  Results  of  18  Tests,  First  Order  Model  of  3  Variables 
Life  predicting  equation: 

y  =  5.9008  -  0.1117x,  -  0.729x„  -  0.013x^ 


Table  10.  Analysis  of  Variance  of  18  Tests 
First  Order  Model  of  8  Variables 
Life  predicting  equation: 

y  =  5.823  +  0.039x,  -  0.904xo  -  0.082x,  +  0.127x. 

12  3  4 

-  0.231X,-  +  0.10 Ox,  -  0.112x,  +  0.015xo 

J  0  /  o 


Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

K-Ratio 

Due  to  Mean 

0.1212 

1 

0.1212 

Due  to  Variance 

8.4720 

1 

8.4720 

Due  to  Zero  Upcrossings 
Due  to  ef  Level 

0.0070 

1 

0.0070 

Upcrossings 

0.1440 

1 

0.1440 

Due  to  Duration  of 
Excursion  Above  Zero 
Duration  of  Excursion 

0.0692 

1 

0.0692 

Above  e,  Level 

Due  to  Band  Width 

0.2291 

1 

0.2291 

0.1988 

1 

0.1988 

Due  to  Average  Amplitude 

Above  ef  Level 

0.0036 

1 

0.0036 

Due  to  Regression 

9.2449 

8 

1.1556 

19.29 

Residual 

0.5394 

9 

0.0599 

Total 

9.7843 

17 

F-ratio  is  greater  than  the  table  value  3.23  with  8  and  9  degrees  of 
freedom  at  95%  significance  level.  So  regression  is  effective  and  the 
model  is  accepted. 


Table  11.  Results  of  18  Tests,  First  Order  Model  of  8  Variables 
Life  predicting  equation: 
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Table  12.  Analysis  of  Variance  of  24  Tests 
First  Order  Model  of  3  Variables 
Life  predicting  equation: 

y  -  5.855  -  0.12lx1  -  0.760x2  -  0.0C6x„ 


Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

F-Ratio 

Due  to  Mean 

0.1172 

1 

0.1172 

Due  to  Variance 

11.5157 

1 

11.5157 

Due  to  Zero  Upcrossings 

0.0021 

1 

0.0021 

Due  to  Regression 

11.6350 

3 

3.8783 

Residuals 

1.4781 

20 

0.0739 

52.48 

Total 

13.1131 

23 

F-ratio  is  greater  than  the  table  value  3.10  with  3  and  20  degrees  of 
freedom  at  95Z  significance  level.  So  regression  is  effective  and 
Che  node!  is  accepted. 
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Table  14.  Analysis  of  Variance  of  24  Tests 
First  Order  Model  of  6  Variables 
Life  predicting  equation: 

y  »  5.826  -  0.037x,  -  0.919x_  -  0.027x„  +  0.078x. 

1  2  _>  4 

-  0.139x5  +  0.053x6  -  0.083x7  +  0.044xg 


Source 

Sum  of 

Degrees  of 

Mean 

F-Ratio 

Squares 

Freedom 

Square 

Due  to  Mean 

0.1172 

1 

0.1172 

Due  to  Variance 

11.5157 

1 

11.5157 

Due  to  Zero  Upcrossings 
Due  to  ef  Level 
Upcrossings 

0.0021 

1 

0.0021 

0.0707 

1 

0.0707 

Due  to  Duration  of 
Excursion  Above  Zero 

Due  to  Duration  of 

0.0554 

1 

0.0554 

Excursion  Above  e. 
Level 

0.2026 

1 

0.2026 

Due  to  Band  Width 

Due  to  Average  Amplitude 

0.1556 

1 

0.1556 

Above  ef  Level 

0.0351 

1 

0.0351 

Due  to  Regression 

12.1545 

8 

1.5193 

23.78 

Residuals 

0.9585 

15 

0.0639 

Total 

13.1131 

23 

F-ratio  is  greater  than  the  table  value  2.64  with  8  and  15  degrees  of 
freedom  at  95%  significance  level.  So  regression  is  effective  and  the 
model  is  accepted. 


Table  15.  Results  of  24  Tests,  First  Order  Model  of  8  Variables.  Life  predicting  equation: 

y  «  5.826  -  0.037x.  -  0.919xo  -  0.027x.  +  0.078x.  -  0.139x,  +  0.053x,  -  0.083x,  +  0.044x, 
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Tabic  16.  Comparison  of  Regression  Sum  of  Squares  for  18  and  24  Tests 


Sum  of 

Squares 

Source 

18  Tests 

24  Tests 

Due  to  Variance 

8.4720 

11.5157 

Due  to  Mean 

Due  to  Duration  of  Excursion 

0.1212 

0.1172 

Above  ef  Level 

0.2291 

0.2026 

Due  to  Band  Width 

0.1988 

0.1556 

Due  to  Regression 

Variance  Alone 

8.4720 

11.5157 

All  4  Significant  Variables 

9.0211 

11.9911 

Total  Sum  of  Squares 

9.7843 

13.1131 

Regression  Sum  of  Squares 
Percent  of  the  Total 


Variance  Alone  86.6  87.8 

All  4  Significant  Variables  92.2  91.4 

Percent  Reduction  in 

Residual  Sum  of  Squares  5.6  3.6 
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Table  17.  Comparison  of  Percent  Deviations  of  Predicted  Lives 
and  Residual  Sum  of  Squares  for  Six  Models 
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